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We present projections for reconstruction of the inflationary potential expected from ESA's up- 
coming Planck Surveyor CMB mission. We focus on the effects that tensor perturbations and the 
presence of non-Gaussianities have on reconstruction efforts in the context of non-canonical inflation 
models. We consider potential constraints for different combinations of detection/null-detection of 
tensors and non-Gaussianities. We perform Markov Chain Monte Carlo and flow analyses on a 
simulated Planck-precision data set to obtain constraints. We find that a failure to detect non- 
Gaussianities precludes a successful inversion of the primordial power spectrum, greatly affecting 
uncertainties, even in the presence of a tensor detection. In the absence of a tensor detection, while 
unable to determine the energy scale of inflation, an observable level of non-Gaussianities provides 
correlations between the errors of the potential parameters, suggesting that constraints might be 
improved for suitable combinations of parameters. Constraints are optimized for a positive detection 
of both tensors and non-Gaussianities. 



I. INTRODUCTION 

A faithful reconstruction of the inflaton potential is 
paramount to understanding the epoch of inflation. Re- 
cent cosmological observations [1-5] have revealed much 
about the possible forms of the potential [6-9]. Future 
CMB missions, such as ESA's Planck Surveyor [10], and 
ground-based polarization experiments such as BICEP 
[11], promise to further hone our understanding of the 
early universe. The Planck Surveyor will provide im- 
proved measurements of the temperature and E-mode 
polarization anisotropics of the CMB out to £ — 3000, 
and will detect the B-mode polarization characteristic 
of gravity waves if the tensor/scalar ratio, r, is greater 
than around 0.05. Planck will also provide the first high- 
quality measurement of any departure from Gaussianity 
exhibited by the CMB temperature fluctuations, with a 
projected sensitivity IJnl] ^ 5 [12, 13]. 

Recent times have also seen exciting progress in in- 
flationary model building, and much is now understood 
about how to implement inflation within string theory 
(see [14-17] for reviews). One particularly exciting real- 
ization is provided by the dynamics of D-branes evolving 
in a higher-dimensional, warped spacetime. A novel as- 
pect of these models is the existence of a speed limit on 
the field space, resulting from causality restrictions on the 
motion of the branes in the bulk spacetime. This allows 
for slow roll to be achieved even in the presence of a steep 
potential. In terms of effective field theory, the speed 
limit is enforced by non-canonical kinetic terms in the 
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DBI action describing the motion of the D-branes. These 
models, termed DBI inflation [18], are further distin- 
guished by the fact that scalar perturbations propagate 
at a sound speed < 1. In the limit ^ 1, these fluc- 
tuations exhibit strong non-Gaussianities [19-23]. Non- 
Gaussianities of this type were first popularized in the 
context of k-inflation [24] , a model- independent construc- 
tion in which higher-derivative kinetic terms drive infla- 
tion in the absence of any potential energy. This phe- 
nomenon is a general characteristic of models with non- 
canonical kinetic terms in the Lagrangian. 

With the advent of these more sophisticated theoret- 
ical constructions, and the increasingly accurate cosmo- 
logical probes planned to study them, more generalized 
methods of inflationary reconstruction have been devel- 
oped. Bean et al. [25] present a formalism for recon- 
structing a general action - both potential and non- 
canonical kinetic energy functions can be linked to ob- 
servations. The new functional degree of freedom pre- 
sented by the non-canonical form of the kinetic energy 
typically gives rise to additional observational signatures. 
For example, the speed of sound, Cs, becomes an ad- 
ditional dynamical degree of freedom in k-inflation and 
DBI, giving rise to non-Gaussian density fluctuations. 
However, Cg also affects the form of the primordial power 
spectra. A successful reconstruction of the inflaton po- 
tential, T^(<^), thus requires a measurement of Cs, since 
otherwise the expressions giving the spectral parame- 
ters cannot be uniquely inverted to obtain V{(j)). An 
accurate reconstruction must also include a detection of 
the tensor power spectrum. Currently, neither of these 
quantities have been measured: present bounds on any 
non-Gaussianities lie within —9 < f^^f^^ < 111 and 
-151 < ff^f < 253 (95% CL) [5], and r < 0.2 (95% CL) 
with WMAP+BAO-KSN [5] and r < 0.35 (95% CL) for 
WMAP-I-ACBAR [26] for power-law spectra, but might 
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be larger if more exotic spectra are considered [27] . 

In this paper, we make reconstruction projections for 
the Planck satellite within the context of this larger in- 
flationary parameter space, namely, one that contains Cg 
as an additional degree of freedom. We seek to deter- 
mine what combination of observations will be required 
to most successfully reconstruct the potential. In the 
best of all worlds, future observations will detect both 
tensors and non-Gaussianities, telling us not only that 
exotic physics is responsible for inflation, but also al- 
lowing for a successful reconstruction program. In this 
paper, we focus particularly on cases in which only one 
of these observables is positively detected. 

In Section 2, we review potential reconstruction for the 
case of non-canonical inflation. In Section 3, we perform 
a Bayesian analysis on a simulated Planck-precision data- 
set for which there is a positive tensor detection, r = 0.1, 
but no detection of non-Gaussianities, I/nl] < 5. This 
corresponds to a universe in which canonical inflation 
took place (with Cg = 1), or, for example, DBI infla- 
tion with Cs € [0.25, 1]. We obtain constraints on V'((/>o), 
V''((/>o), and V"{(j)o) for both possibilities and find that 
the error bars arc significantly larger if Cg is taken as a 
free parameter within the above range. Our inability to 
resolve Cg prohibits us from using the shape of the power 
spectra alone to reconstruct V{(l)). In Section 4, we inves- 
tigate this problem using infiationary fiow methods. The 
flow method is a natural extension of the MCMC analy- 
sis, allowing higher dimensional parameter spaces to be 
explored with relative ease. We reproduce the results 
obtained in Section 3, and then investigate the effect of 
allowing the speed of sound to vary over the course of 
inflation. We find that constraints on Vq loosen consid- 
erably as this new degree of freedom is introduced. Next, 
we use flow methods to investigate the observational pos- 
sibility of a positive detection of /nl and a null detection 
of r. We find that constraints on the individual potential 
parameters are worse in this case than in the case where 
tensors are detected, but that the errors are highly cor- 
related, and tight bounds can be placed on combinations 
of these parameters. Finally, we consider the case for 
which both r and Jnl are measured. This is clearly the 
best case scenario for near future experiments, since it 
will not only place bounds on the potential parameters 
from above and below, but it will also be a smoking gun 
for exotic physics. In Section 5 we present conclusions. 



II. INFLATION MODELS WITH AN 
ARBITRARY SPEED OF SOUND 

A. Flow parameters and reconstruction 

In this section we consider the case of a scalar field 
with a general Lagrangian C{X, cj)), where (f) is the field 

that is responsible for inflation {the inflaton) and X = 
(1/2)5'"' If the inflaton field is coupled to Ein- 



stein gravity, the total action S can be written [28] 
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where R is the curvature scalar and we work in units of 
the reduced Planck mass, Mpi = mpi/\/8n. We can also 
describe the inflaton field ^ as a, perfect fluid by express- 
ing the energy- momentun tensor T^^ in the following way 
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where the pressure P is identified with the Lagrangian 
jC{X, (j)), p is the energy density of the field, 
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is the fluid 4-velocity. Assuming that the four- 
dimensional metric is of the Friedmann-Robertson- 
Walker (FRW) form 

Qf,, = dmg{-l,a\t),a^(t),a^(t)) , (5) 

the Friedmann and the acceleration equations are [25] 
2 
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The speed at which fluctuations propagate relative to 
the homogeneous background is determined by the sound 
speed of the fluid. 



dP 



P. 



X 



1 



(7) 



- dp P,x+2XP,xx 7' 
Infiation models driven by canonical scalar fields satisfy 

P,xx = 0, (8) 

and so fluctuations travel at the speed of light 

= 1. (9) 

However, in general non-canonical models, fluctuations 
can propagate at speeds much lower, resulting in non- 
Gaussian adiabatic fluctuations [19-22]. In terms of the 
comoving curvature perturbation, the existence of non- 
Gaussianities is reflected in a non-trivial three-point func- 
tion, (CkiCk2Ck3)- In non-canonical models the dominant 
source of non-Gaussianity is due to configurations for 
which the wavenumbers ki = k2 = k^, forming equi- 
lateral triangles in Fourier space. To lowest order in slow 
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roll, the magnitude of this non-Guassianity is given by 

the parameter [20], 
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with 
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In the context of more general actions such as Eq. (1), 
the process of reconstructing the physics of inflation re- 
quires an understanding of both the potential energy of 
the field, V {(/)), and the kinetic term. Such a formalism 
was recently developed in Ref. [25] in terms of different 
hierarchies of flow parameters. Such parameters natu- 
rally arise in the Taylor expansions of the functions H, 
7, and C,x, which together completely specify the in- 
flationary dynamics. In this paper, we fix the gauge by 
setting [29] 
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which leads to the following hierarchies of flow parame- 
ters [30], 
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where £ = 2, • • • , oo is an integer index. For the remain- 
der of this paper, we will refer to the parameters ^X{(f)) 
and ^a{<j)) as the £-th components of the H-tower and 
^ -tower respectively. In the case where 7 = 1, the 7- 
tower vanishes and the above hierarchy reduces to the 
flow hierarchy of a canonical scalar field [31, 32]. 

We next Taylor expand the Hubble parameter about 
<^o = 0, 



H{cP) = Ho 



and the inverse speed of sound 
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where the expansions are truncated at orders -|- 1 and 

Af +1, respectively. 
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for £>M + 2, 



= for £>J\f + 2. 
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In fact, this truncation is equivalent to simply setting 
-^+1A = ^~^^a = in the flow hierarchy, for reasons 
that will become clear in Section IV. 

The coefficients A( and B( are given by [30] 
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While this particular parameterization restricts H{(j)) 
and 7(0) to be polynomials, these functions can be made 
more general by truncating Eq. (13) at arbitrarily high 
order. The values of the flow parameters at some 4> = 4>o 
fully define the functions H{(f)) and completely de- 

termining the theory in the gauge £,x = cj^. 



B. DBI Inflation 

In what follows, we will study DBI inflation as a pro- 
totype non-canonical model. The only model-dependent 
results in this paper are those relating non-Gaussianities 
to the sound speed, c^. However, as we will discuss fur- 
ther, the qualitative nature of our results should apply 
generically to non-canonical models. DBI inflation refers 
to models arising from the action of D-branes in flux 
compactifications of type-IIB string theory [18, 33, 34]. 
Typically, the inflaton field parameterizes the separation 
of a probe D3-brane and a stack of D3-branes in a warped 
throat geometry. This geometry is defined by the line el- 
ement [35], 



ds 
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(19) 

where the internal space is a cone over the five-manifold 
X5. The inflaton field is proportional to the throat co- 
ordinate /3 by (/) = \/T^p, where T3 is the tension on the 
D3-brane. The Lagrangian for the inflaton is then of the 
form [34] 



= -/-i(0)0T/(0)rt^ + r\^) - V{^), 

(20) 
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where V{(l)) is the scalar field potential and the function 
/((/)) is related to the warp factor h{(p) by 



m 

The 7-factor defined by 
7 = 



1 



(21) 



(22) 



enforces a speed limit on the moduli space, arising from 
causality constraints on the motion of the probe brane 
in the throat. This phenomenon mitigates the //-problem 
associated with potentials derived from supergravity, al- 
lowing for otherwise too steep potentials to successfully 
drive inflation. 

The equations of motion are easily obtained from our 
previous results. The Friedmann equation (6) becomes 
[36], 



\<t>) - V{ct^) 



which gives 



7(0) 

Equation (22) then becomes, 
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7 = ^Jl + 4M^,m[H'{cP)f. (25) 
These results lead to the Hamilton-Jacobi equation, 



V{(l)) = 3M|i/f2((/)) - 4M^i 
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which can be used to reconstruct the scalar potential 
given the functions i?((/>) and 7(<^). The warp factor of 
the geometry can also be reconstructed from these two 
fmictions via Eq. (2-5). The Taylor expansions Eq. (14) 
and (15) can then be used to obtain exact expressions 
for V{<l)o) and its derivatives. In this work, we will only 
require expressions for the first two derivatives of V{(j)). 
After some tedious but straightforward algebra we ob- 
tain, 

ViM = M^iH^ {3 - 2eT) , 



Mpiij2 V2e7 (3 - 277r + sT^) 
+2s{2r]-sr)T^ + epT^] , 



where 
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In the next section, we utilize these results to obtain ob- 
servational constrains on the inflaton potential from a 
synthetic data-set. 



III. PROJECTIONS FOR PLANCK: MCMC 
ANALYSIS 

In this section, we constrain the inflaton potential 
in the hypothetical case that Planck detects tensors, 
with r = 0.1, but fails to positively detect any non- 
Gaussianities, implying |/jvl| ^ 5. Such an observation 
would certainly be consistent with canonical inflation, 
for which Cg = 1. However, it would also be consistent 
with a host of more exotic, non-canonical models, such 
as k-inflation and DBI. The speed of sound varies in such 
models and gives rise to detectable non-Gaussianities if 
Cs <C 1. From Eq. (10), we have for the case of DBI 
inflation, 



f - ^5 / 1 
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However, even for non-detectable values of /nl, we see 
that the sound speed can vary considerably, 0.25 ^ Cg < 
1. In non-canonical models, there is another way to mea- 
sure Cg, however, independent of the strength of non- 
Gaussianities. These models predict a modified consis- 
tency relation [28], 



(30) 



where rit is the tensor spectral index. The downside is 
that this relation will be very difficult to reliably measure 

any time in the near future. Even a direct detection 
of primordial gravity waves from proposed BBO [37] or 
DECIGO [38] will only measure this relation to within 

an accuracy of 50% [39]. The uncertainty in the value 
of Cs thus remains. How might this affect our ability to 
reconstruct V(0)? 

In models with a variable speed of sound, the scalar 
and tensor amplitudes are given during slow roll by. 



1 



kcs—aH 



and 



2 H 



PI 



(31) 



(32) 



k=aH 



where we note that scalar perturbations freeze out upon 
exiting the sound horizon, whereas tensor perturbations 
do so upon exiting the Hubble radius. Often in the lit- 
erature, the tensor/scalar ratio is written r = 16cse, ob- 
tained from Eqs. (31) and (32) by neglecting this differ- 
ence in freeze out times (assuming that H is constant). 
This is not strictly correct. While this difference enters at 
2"''-order in slow roll, it can be important if Cg is small^ . 



^ The limited accuracy of this expression was recently noted in Ref. 
[29] , where numerical spectrum calculations reveal discrepancies 
of up to 50%. 
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The evolution of H is determined by 



(33) 



where N is defined as the number of e-folds before the 
end of inflation. Prom the horizon crossing conditions we 
form the ratio, 



c, at 



(34) 



where the subscripts t and s indicate that these values are 
measured when the tensor and scalar perturbations exit 
their respective horizons. If we suppose that e is constant 
between these crossing times (an excellent assumption 
during slow roll), this ratio becomes 



(35) 



with the result that the change in the number of e-folds 

between the time that a certain scalar fc-modc exits the 
sound horizon and the time the corresponding tensor 
mode exits the event horizon is given by. 



iV. -N,= 



The tensor /scalar ratio is then' 
Pn{k) - '^'"^ U. 
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(36) 
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The scalar spectral index of the power spectrum is also 
affected by the variable speed of sound. At lowest-order 
in slow roll, 



1 - 4e + 2r/ - 2s, 



(38) 



where rj and s are defined in Eq. (13). The tensor spec- 
tral index is written 



nt = -2e, 



(39) 



which leads to the modified consistency relation at lowest 
order, r = —ScsTit. 

We are now in a position to determine what effect a 
non-detection of Jnl has on reconstruction efforts. Prom 
our results in the previous section Eq. (27) we can link 
V{<j)) to the observables P-jz, r, Ug, and /nl- In this 
section we work to lowest-order in slow roll and assume 
that Cs is constant. Wc neglect any running of the spec- 
tra and thus work to reconstruct V{(j)o) and its first two 
derivatives. In terms of observables these are written. 
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where the observables are to be measured at the scale 
corresponding to (po- We have used 7 rather than the 
observable /jvl in these expressions for simplicity. It is 
now clear how any uncertainty in 7 will translate into un- 
certainties in the potential. The height of the potential, 
V{(j)o), will be least affected, however, V'{(j)o) and V"{(po) 
will be more strongly so because they are directly pro- 
portional to 7. 

In order to rigorously analyze the effects of an un- 
resolved speed of sound, we constrain l^(0o) and its 
first two derivatives using Markov Chain Monte Carlo 
(MCMC) and a simulated, Planck-precision CMB data- 
set. We use the publicly available CosmoMC'^ [41] package 
to explore the inflationary parameter space. We con- 
strain V{(j)) arising from non-canonical models by allow- 
ing Cs to vary in the analysis, and compare this case to 
constraints obtained on V{(j)) arising from canonical in- 
flation, obtained by fixing Cg = 1. 

MCMC techniques [42-45] sample the likelihood sur- 
face of model parameters, C{d\6), where d represents the 
n-dimensional data and the n-dimensional parameter 
vector. The likelihood function relates any prior knowl- 
edge of the parameter values, n{6), to the posterior prob- 
ability distribution via Bayes' theorem. 
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The posterior probability distribution function (PDF) of 
a single parameter Oi can then be obtained by marginal- 
izing p{6\d) over the remaining parameters, 



p{ei\d) = j p{e\d)d9i ■ ■ ■ de^-ld9i+^ ■ ■ ■ dOn- 



(42) 



In this study, we only constrain the inflationary param- 
eters {V{(j)), V'{<j)), V"{(j)), and Cg), fixing the late-time 
parameters at the values: baryon and CDM densities, 
Qbh"^ = 0.022 and Qch"^ = 0.104, the angular diameter 
distance at decoupling, 0g = 1.04, and the optical depth 
to reionization, r = 0.093. Degeneracies amongst the 
late-time parameters and the power spectrum do exist, 
in particular, n^ is still moderately degenerate with flj,h^. 
However, our goal in this study is to compare the con- 
straints obtained on V{<j)) when Cg is allowed to vary and 
when it is held constant at = 1. In both of these cases 
the spectrum is completely described by the parameters 



^ See Ref. [40] for an alternative derivation. 



^ http;//cosmoIogist.info/cosmomc/ 
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FIG. 1: Marginalized 1- and 2-D PDF's of the potential parameters, Vo = V{(l>o), Vi = V'{(l>o), and V2 = ^^"(9^0) in units 
of Mpi. The small red contours correspond to 1- and 2-cr errors in canonical inflation, in which Cs is held fixed at 1. The 
substantially larger black contours correspond to non-canonical infiation in which Cs is allowed to vary in the chains between 
0.25 and 1. 



Ptz, r, ris, and nt, and so the introduction of Cs as a free 
parameter will not introduce new degeneracies. 

Our Planck-precision simulation comprises tempera- 
ture, E- and B-mode polarization data out to £ = 2500 
with a sky coverage of 65%. For simplicity, we consider a 
single-channel experiment in which the full width at half 
maximum of the Gaussian beam, ^twhm = 7', and the root 
mean square of the pixel noise, = ap/2 ~ 30.25/iK^. 
As a fiducial model we choose a power-law spectrum with 
Pn{k = 0.002/1 Mpc"^) = 2.3 x 10"^, = 0.97, and 
r — 0.1. We assume the inflationary consistency rela- 
tion for the tensor spectral index, nt — —"fr/8. We im- 
pose the prior Cg G [0.25, 1], consistent with undetectable 
non-Gaussianities. We assume purely adiabatic initial 
perturbations and spatial flatness, and adopt a top-hat 
prior on the age of the universe: 10 < to < 20 Gyr. We 
utilize the Metropolis-Hastings sampling algorithm and 
measure convergence across 8 chains with the Gelman- 
Rubin R statistic. 

We present our results in Figure 1 and Table 1. 
Both best-fit reconstructions fit the data equally well: 
-21n£ = 7931.29 for the canonical and -2lnC = 7931.16 
for the non-canonical reconstruction. The two models 
are thus statistically degenerate. The red contours in the 
figure correspond to 1- and 2-a errors for canonical infla- 
tion, for which Cs is held flxed at 1. The black contours 



correspond to non-canonical inflation in which Cs is al- 
lowed to vary in the chains between 0.25 and 1. Indeed, 
the constraints on V{(j)) worsen significantly if Cs is in- 
cluded in the reconstruction. This is a significant imped- 
iment to the program of potential reconstruction: if we 
allow for the possibility of an arbitrary speed of sound, 
reconstruction of the potential is severely degraded in 
the absence of a measurement of the amplitude of non- 
Gaussianities. Table I indicates that the marginalized 
errors on Vi and V2 increase by an order of magnitude in 
going from canonical to non-canonical inflation. While 
we have used Eq. (29) specific to DBI infiation, the fact 
that a failure to resolve non-Gaussianities precludes a 
unique inversion of the power spectrum is a general re- 
sult. We therefore expect this problem to persist in gen- 
eral non-canonical models, however, the degree of the 
effect is model dependent. 

It is interesting to note how well constrained the canon- 
ical model is. In particular, the absolute errors on V" are 
smaller than those obtained on V' . The dominant term 
in the error on V" is proportional to 3Tr^drP'ji{ns — l)/2, 
whereas the dominant term in the error on V' is SV' ~ 
Pn{r + i5r)3/2 - V'. For Planck -precision errors on r, 
the two terms in SV are of the same order, and so 
0{5V') ^ 0{V'). For the values of our fiducial model. 
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Model 


Vo X lO'^Mp/ 


Vi X lO^Mpi^* 


V2 X 10"Afpi^ 


canonical 


n 7+1.9 
' -0.7 




^•^-1.8 


non-canonical 


n g+1.9 
■^■"-0.9 


^"-12 


01+80 
^^-90 



TABLE I: Marginalized 2-a errors on the potential parameters 
of each model. 



(r + Jr)"^/^ « 5r, and so roughly we have 

5V" « ^Ti'^in, - 1)5V' ~ O(10-i)<5F'. (43) 

In Figure 2 we plot the reconstructed best fit poten- 
tials for each model. As indicated by the error contours 
and by this figure, the speed of sound has the effect of 
rendering steep potentials viable. This is precisely the 
novel aspect that initially popularized DBI inflation, but 
as we see here, is a bane to reconstruction efforts if /at^ 
goes undetected. However, the detection of tensor modes 
still reliably determines the energy scale of inflation. For 
r = O.I, the uncertainty introduced by Cs is only around 
1%. The energy scale reconstructed from the observables 
is measured at = in Figure 2. 



IV. PROJECTIONS FOR PLANCK: FLOW 
ANALYSIS 

While the MCMC method utilized in the last section 
is a rigorous Bayesian analysis of the parameter space, 
it performs poorly if too many free parameters are in- 
cluded in the fit. This is especially true if some of these 
parameters are degenerate or poorly constrained by the 
data. It is therefore not feasible to go much beyond the 
lowest-order analysis (cs ~ const.) conducted in the pre- 
vious section. In this section, we use the flow formalism 
[32, 46] as a method for reconstructing inflationary po- 
tentials compatible with future observations. While not 
statistically rigorous, the flow method is well suited to 
studying the inflationary parameter space to arbitrarily 
high order in slow roll. It is also capable of reconstruct- 
ing models with a varying sound speed. In this section, 
we utilize the flow formalism in order to study parameter 
constraints in the event that a tensor signal is detected 
in the future, but non-Gaussianities are not. This is the 
same scenario considered in the last section, but here we 
examine the effects of allowing a time-varying speed of 
sound. We conclude this section by considering the case 
of a positive detection of non-Gaussianities. 

The flow formalism is based on the property of the flow 
parameters defined in Eq. (13), that their evolution in 
terms of the number of e-folds N can be described by a 
set of first order differential equations. Taking successive 
derivatives of the flow parameters with respect to iV, we 
obtain the flow equations 

I dH 
^ " HdN' 
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FIG. 2: Potentials reconstructed from the best-fit parameters 
of Table 1. Note the increased steepness of the non-canonical 
potential (black dashed) relative to the canonical potential 
(red solid). 



dt 

dN 
drj 

dN 



fx 

dN 



ds 

dN 
dp 

dN 



e(2r?-2e-s), 
-?]{e + s) + ^X, 

I d7 
7d/V' 

—s {2s + e ~ rj) + ep, 



= -2ps+ a, 



— = -^a[(f+I). + (€-I)(e-,?)]+^+V(44) 

where the derivative of a flow parameter is related to 
terms of higher-order in the flow hierarchy. In practice, 
the flow equations are truncated at some flnite order M 
and Af by requiring that = = for all £ > + I 
and Af +1, respectively, at the initial time. By Eq. (44), 
these remain zero for all time with the result that the tra- 
jectory so obtained is an exact solution to the inflationary 
equations of motion. 

This truncated system of differential equa- 
tions can be solved numerically by specifying the 
set of initial conditions for the flow parameters 



[H, e, . . . A; 7, s, . . . ,-'^a\ at some time Nq. The 
solution will then be a particular inflationary path 
[HiN), eiN), . . . ,^ X{Ny,^{N),s{N), . . . ,^a{N)] in 
the (A4 -(- A/')-dimensional parameter space. In fact, 
the expansions Eqs. (14) and (15) yield exact analytic 
solutions of the flow equations, where the expansions 
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can be taken about any point along the inflationary 
trajectory [47]. In this paper, we truncate the flow 
hierarchy to second order in the iJ-tower, and to first 
order in the 7-tower so that 



(45) 



for all £ > 2 and £' > 1. The resulting solu- 
tion to the flow equations will then be a specific path 
[e{N),r]{N),^{N),s{N),p{N)] in the 5-dimensional pa- 
rameter space. The inclusion of ,s and p allows us to 
generate models with a varying speed of sound. We solve 
the flow equations from the time the quadrupole leaves 
the horizon at iVquad up until the smallest scale modes 
detectable in the CMB are generated, a period typically 
spanning 10 e-folds. We restrict this analysis to models 
that do not violate causality, < 1. We then calculate 
the observables from the values of the flow parameters at 

N = A/-quad, 



ris = 



1 - 4e + 2r/ - 2s - 2(1 + C)e'^ - (3 + C)s'^ 
-^(3 - 5C)er? - ^(11 + 3C)es + (1 + C)r?s 

+^(l + C)6p+i(3-C)(2A), 
16c.e l + l(e-,?)(C-3) + |(C + l)| , 



1 



where 



1-e-sJ dN' 

35~ / 1 



108 



1 



C = 4(ln2 + 7) - 5, 



(46) 



(47) 



and 7 ~ 0.577 is the Euler-Mascheroni constant, which 
should not be confused with the inverse sound speed Eq. 
(22). 



do not have analytic expressions'* for these higher-order 
terms, we can not be sure that they are under control. 
By working only up to order and by ensuring that 
the running is small, we can be certain that all higher- 
order terms will likewise be small, since they are at least 
of order Most viable models at this order naturally 
have small running anyway, since large negative running 
at order typically leads to an insufficient amount of 
inflation [48]. 

We now use this methodology to stochastically gener- 
ate models of inflation consistent with the constraints Eq. 
(48). We are interested in models which fail to generate 
observable non-Gaussianities, corresponding to a speed 
of sound in the range G [0.25, 1] (c.f. Eq. (29)). We 
first seek to reconstruct models with a constant speed of 
sound in order to make contact with the results of section 
3, and so the 7-tower is simply replaced by Cs = const. 
After obtaining results for constant speed of sound, we 
perform the analysis for a variable speed of sound by 
allowing s and p to vary. 

We present our results in Figure 3. We plot the val- 
ues of Vo, Vi = Vq, and V2 = Vq for canonical inflation, 
where Cg = 1 (red squares), DEI inflation with constant 
speed of sound (black circles), and DEI inflation with 
varying speed of sound (green crosses) using Eqs. (27). It 
can easily be seen that, even though we can not make any 
rigorous statistical arguments using the flow formalism, 
we successfully reproduce the results found in the previ- 
ous section using MCMC techniques. The most striking 
consequence of a variable speed of sound is that a large 
region of parameter space for which Vq < opens up. 
To understand this, note from Eq. (27) that the sign of 
V" to lowest-order is determined by the sign and mag- 
nitude of rj. In the case of constant sound speed, the 
value of r] is tightly constrained by the spectral index, 
Ug = 1 — 4:6 + 2-q. However, if Cg is allowed to vary, then 
the spectral index becomes 



l-Ae + 2'q- 2s, 



(49) 



A. Detection of r and no detection of /ivr 

We are interested in collecting models compatible with 
the simulated data of Section III, 



r 
a 



0.97 ±0.0036, 
0.1 ±0.05, 
0.0 ±0.005, 
2.3 ±0.102, 



(48) 



to lowest-order. For a given value of n,,, rj can be taken 
more negative by suitably adjusting s. From Eq. (27), 
this allows a wide range of parameter space for which 
Vq < to be brought into agreement with the data. 
Physically, the steepening of the potential that results 
from tuning r] is mitigated by the increased warping of 
the geometry, which is related to the value of the flow 
parameter s. 



where the errors are l-u projections for Planck, quoted 
at fc = 0.002 /iMpc~^. The observables calculated at the 
quadrupole by the fiow method are easily extrapolated to 
k = 0.002 /iMpc~^ because the spectra are power-laws. 
In order to efficiently generate models with power-law 
spectra, we truncate the ff-tower at The inclusion 
of higher-order terms not only aff'ects the running, but 
may even make a non-negligible contribution to higher- 
order fc-dependencies in the power spectrum. Since we 



B. Detection of /ati. 

In this section, we are interested in models which are 
characterized by detectable non-Gaussianity for both ob- 



^ This issue is avoided by calculating the power spectra numeri- 
cally, as done in [29]. 
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FIG. 3: Results of the flow analysis for the case of detection of r, with r — 0.1 ± 0.05 but no /jv_t {/nl < 5). Red squares are 
canonical models, black circles DBI inflation with constant speed of sound, and green crosses DBI models with variable speed 
of sound. The case of variable speed of sound opens up the parameter space for which Vq" < 0. 



servational outcomes of gravity waves. We start by con- 
sidering the case in which future data fails to detect any 
gravity waves, but does successfully detect evidence for 
non-Gaussian fluctuations. While it is well known that 
a failure to detect tensors precludes a lower bound on 
the energy scale of inflation, it might be hoped that a 
positive detection of non-Gaussianities yields improved 
constraints relative to the worst case scenario in which 
neither is detected. The simulated data used thus far in 
the study does not apply to this case. However, taking 
the scalar perturbations to be the same, we collect mod- 
els that satisfy Eq. (48) for Ug, Pn, and a. In order for 
Planck to fail to detect tensors, we suppose that the null 
hypothesis (r = 0) is verified for r < 0.05. 

We show the results of the flow analysis in Figure 4. We 
consider detection of non-Gaussianities at two different 
levels: /nl = 10 and Jnl — 40, with l-tr errors AfffL = 
±5. By failing to detect r, we only succeed in placing 
an upper bound on the values of Vq and Vq. However, 
a determination of 7 tightly correlates Vq with both Vq 
and Vq. From Eq. (41) we have that 



V 



Pn{k) 



1/3 



y/2/3 



and 



V 



Pn{k) 

1 




AV" 



1 



(50) 



(51) 



The uncertainty in r only allows us to constrain the above 
relationships between the potential parameters. These 
relations give the distinctive band structure seen in Fig- 
ure 4. The effect of the detection of non-Gaussianities on 
these constraints is two-fold: the value of 7 controls the 



'slope' of the bands, and, as we show next, the relative 
error between V-V and V-V" . 

The contribution of the error A/nl to V and V" de- 
pends on the absolute value of /nl- Since this is the 
dominant source of error, a determination of /nl can 
thus improve constraints. To see this, consider the case 
of V . Including only the contribution of /jvl to the error 
on V' , we have 

SV CX (V/iVL + AfNL - VfNL - A/nl) SV, (52) 

and expanding the square roots for small A/jvl relative 
to Jnl gives 



Sv'o^^SV 



V7 



NL 



(53) 



Thus, a determination of /nl correlates the errors in V 
and V, and similarly for V and V" (c.f. Figure 4). Addi- 
tionally, the larger the value of Jnl, the smaller this rel- 
ative error between the potential coefficients. Therefore, 
while a detection of r allows us to constrain the absolute 
scale of the inflationary potential and its derivatives, a 
detection of /nl gives us information about the relative 
values of these parameters, as well as being a smoking 
gun for exotic physics. 

Finally, we consider the case where both r and f^^ 
are detected. We collect models that satisfy Eq. (48), 
and produce detectable non-Gaussianities at the level of 
Jnl = 10 with 1-a errors Af^L = ±5. We present the 
results of the flow analysis in Figure 5, where it can be 
seen that we again obtain the band-shaped distributions 
for the potential parameters. Following the reasoning 
of the previous sections, we see that the detection of r 
places bounds on the potential parametes, where the de- 
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FIG. 4: Results of the flow analysis for the case of detection of /nl but no r (r < 0.05). The top figure corresponds to models 
with /ml ~ 10 ± 5, and the bottom one to models with /nl = 40 ± 5. Black circles are DBI models with constant speed of 
sound, and green crosses DBI models with variable speed of sound. 



tection of non-Gaussianities determines the relative error 
amongst them. 



V. CONCLUSIONS 

We have considered potential reconstruction forecasts 
for the upcoming Planck Surveyor CMB mission in the 
context of non-canonical inflation. We have focused on 
the following observational outcomes: a detection of a 
tensor signal and a null detection of non-Gaussianities, 
a null detection of tensors and a positive detection of 
non-Gaussianities, and a positive detection of both ten- 
sors and non-Gaussianities. To explore the first case, we 
perform a Markov Chain Monte Carlo analysis on a sim- 
ulated Planck-precision data-set. We obtain constraints 
on the potential parameters Vq, 1q', and Vq , and the 
sound speed Cg. A null detection of non-Gaussianities 
at Planck resolution corresponds to I/atl] ^ 5, with the 



result that c<, is only constrained to lie within the range 
[0.25, 1] for the case of DBI inflation. Since Cg also affects 
the form of the power spectrum, a failure to constrain it 
precludes us from inverting our spectrum observation to 
obtain the inflaton potential. The errors on Vq and Vq 
increase by an order of magnitude relative to constraints 
obtained on canonical inflation. 

We next used the flow formalism to investigate the ef- 
fect of allowing the speed of sound to vary over the course 
of inflation. This requires the introduction of terms gov- 
erning its evolution. Unlike MCMC methods, the flow 
formalism is well suited to studying reconstructions to 
high-order. The strongest effect of a variable speed of 
sound is a much increased error on the second derivative 
of V{(j)). A variable speed is generic to models such as 
DBI inflation. The worsening of the constraint on Vq 
results from the fact that there is additional freedom in 
tuning the geometry in order to get the same observ- 
ables. (We note that this has been a very conservative 



11 




1 ' I ' I ' r 




I ■ I ■ I ■ I ■ I . I 



6e-09 

5e-09 

4e-09 

3e-09 ^ 
> 

2e-09 
1e-09 





1e-09 2e-09 3e-09 4e-09 5e-09 -1e-09 1e-09 2e-09 3e-09 4e-09 5e-09 
Vi((^) V2((^) 



FIG. 5: Results of the flow analysis for the case of detection of r, with r = 0.1 ± 0.05 and /jvl, with /nl = 10 ± 5. Black 
circles are DBI models with constant speed of sound, and green crosses DEI models with variable speed of sound. 



analysis as far as our constraints on f^i^ are concerned. 
If one considers /^^'\ the non-Gaussianities associated 
with equilateral triangles in /c-space, the Planck detection 
threshold could be as large as f^f < 60 [49].) 

We also used the flow approach to investigate the ob- 
servational outcome of a detection of different levels of 
non-Gaussianity. We first consider the case of a null 
measurement of r. Because neither Vq or Vq are bounded 
from below, and since a larger level of non-Gaussianity 
is generated by larger values of VJ, we find that the er- 
rors on the individual parameters Vq, Vq, and Vq grow 
relative to the case where neither are detected. However, 
these errors are strongly correlated, and it is nonethe- 
less possible to tightly constrain combinations of these 
parameters. The parameters r and /nl therefore play 
complementary roles in potential reconstruction. Finally, 
we consider the case where both r and Jnl sue detected. 
The complementary role that r and l play in potential 
reconstruction is again clear, where gravity waves place 
bounds on Vq and Vq from above and below, and non- 
Gaussianities determine the width of the bands in the 
V-V and V-V" planes. 

The approach taken in this paper has been phenomeno- 
logical. We have treated the potential and Cs as freely 
tunable parameters. We have purposefully kept the study 



general to encompass the larger class of fc-inflation and 
other non-canonical models based on effective field the- 
ory. However, it is very likely that many of the potentials 
and sound speeds sampled in this analysis do not repre- 
sent realistic compactification scenarios in string theory, 
and might not correspond to realizations of DBI infla- 
tion. What we have shown in this study is that, if the 
observational outcomes considered here are reflected in 
the actual data, then our ability to reconstruct the po- 
tential is signiflcantly weakened. While we ultimately 
hope that the best-case scenario will be realized with a 
detection of both tensors and non-Gaussianities, in the 
meantime we can only wait with guarded anticipation as 
to what Planck will reveal about the universe. 



Acknowledgments 

We acknowledge the use of the UB Physics Graduate 
Computing Facility and thank Mark Kimball for assis- 
tance. This research is supported in part by the National 
Science Foundation under grants NSF-PHY-0456777 and 
NSF-PHY-0757693. This work was supported by World 
Premier International Research Center Initiative (WPI 
Initiative), MEXT, Japan. 



[1] G. Hinshaw et al. [WMAP Collaboration], 

arXiv:0803.0732 [astro-ph]. 
[2] B. Gold et al. [WMAP Collaboration], arXiv;0803.0715 

[astro-ph] . 

[3] M. R. Nolta et al. [WMAP CoUaboration] , 

arXiv:0803.0593 [astro-ph]. 

[4] J. Dunkley et al. [WMAP Collaboration], 



arXiv:0803.0586 [astro-ph]. 
[5] E. Komatsu et al. [WMAP Collaboration], 

arXiv:0803.0547 [astro-ph]. 
[6] W. H. Kinney, E. W. Kolb, A. Melchiorri and A. Ri- 

otto, Phys. Rev. D 74, 023502 (2006) [arXiv:astro- 

ph/0605338]. 

[7] J. Lesgourgues and W. Valkenburg, Phys. Rev. D 75, 



12 



123519 (2007) [arXiv:astro-ph/0703625]. 
H. V. Pciris and R. Easthcr, JCAP 0807, 024 (2008) 
[arXiv:0805.2154 [astro-ph]]. [31 
L. Lorenz, J. Martin and C. Ringeval, Phys. Rev. D 78, 
063543 (2008) [arXiv:0807.2414 [astro-ph]]. [32 
See http : / /www . rssd . esa . int/ index . php?pro j ect=PLANCK. 
K. W. Yoon et al., arXiv:astro-ph/0606278. [33 
E. Komatsu, arXiv:astro-ph/0206039. 
P. Creminelli, JCAP 0310, 003 (2003) [arXiv:astro- [34 
ph/0306122]. 

J. M. Cline, arXiv:hep-th/0612129. 

R. Kallosh, Lect. Notes Phys. 738, 119 (2008) [35 

[arXiv:hep-th/0702059]. 

C. P. Burgess, PoS P2GC, 008 (2006) [Class. Quant. [36 
Grav. 24, S795 (2007)] [arXiv:0708.2865 [hop-th]]. 
L. McAUister and E. Silverstein, Gen. Rel. Grav. 40, 565 [37 
(2008) [arXiv:0710.2951 [hep-th]]. [38^ 
E. Silverstein and D. Tong, Phys. Rev. D 70, 103505 
(2004) [arXiv:hep-th/0310221]. [39 
M. Alishahiha, E. Silverstein and D. Tong, Phys. Rev. D 
70, 123505 (2004) [arXiv:hep-th/0404084]. [40 
X. Chen, M. X. Huang, S. Kachru and G. Shiu, JCAP 
0701, 002 (2007) [arXiv:hep-th/0605045]. [41 
M. SpaUnski, Phys. Lett. B 650, 313 (2007) [arXiv:hep- 
th/0703248]. [42 
R. Bean, X. Chen, H. V. Peiris and J. Xu, [43 
arXiv:0710.1812 [hep-th]. 

M. LoVerde, A. Miller, S. Shandera and L. Verde, [44 
arXiv:0711.4126 [astro-ph]. 

C. Armendariz-Picon, T. Damour and V. F. Mukhanov, [45 
Phys. Lett. B 458, 209 (1999) [arXiv:hep-th/9904075] . 
R. Bean, D. J. H. Chung and G. Geshnizjani, 
arXiv:0801.0742 [astro-ph]. [46 
W. H. Kinucy, E. W. Kolb, A. Melchiorri and A. Riotto, 
arXiv:0805.2966 [astro-ph]. [47 
B. A. Powell and W. H. Kinney, JCAP 0708, 006 (2007) 
[arXiv:0706.1982 [astro-ph]]. [48; 
J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 219 
(1999) [arXiv:hcp-th/9904176]. [49; 
N. Agarwal and R. Bean, arXiv:0809.2798 [astro-ph]. 
H. V. Peiris, D. Baumann, B. Friedman and A. Cooray, 



Phys. Rev. D 76, 103517 (2007) [arXiv:0706.1240 [astro- 
ph]']. 

A. R. Liddlo, P. Parsons and J. D. Barrow, Phys. Rev. D 
50, 7222 (1994) [arXiv:astro-ph/9408015]. 
W. H. Kinney, Phys. Rev. D 66, 083508 (2002) 
[arXiv:astro-ph/0206032] . 

S. Kachru, R. Kallosh, A. Linde and S. P. Trivedi, Phys. 
Rev. D 68, 046005 (2003) [arXiv:hep-th/0301240]. 
S. Kachru, R. Kallosh, A. Linde, J. M. Maldacena, 
L. P. McAllister and S. P. Trivedi, JCAP 0310, 013 
(2003) [arXiv:hep-th/0308055]. 

M. R. Douglas and S. Kachru, Rev. Mod. Phys. 79, 733 
(2007) [arXiv:hep-th/0610102]. 

M. Spalinski, JCAP 0704, 018 (2007) [arXiv:hep- 
th/0702118]. 

See http:/ /universe. nasa.gov/new/program/bbo. html 

N. Seto, S. Kawamura and T. Nakamura, Phys. Rev. 

Lett. 87, 221103 (2001) [arXiviastro-ph/OlOSOllJ. 

T. L. Smith, H. V. Peiris and A. Cooray, Phys. Rev. D 

73, 123503 (2006) [arXiv:astro-ph/0602137]. 

L. Lorenz, J. Martin and C. Ringeval, Phys. Rev. D 78, 

083513 (2008) [arXiv:0807.3037 [astro-ph]]. 

A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002) 

[arXiv:astro-ph/0205436] . 

N. Christensen and R. Meyer, arXiv:astro-ph/0006401. 
N. Christensen, R. Meyer, L. Knox and B. Luey, Class. 
Quant. Grav. 18, 2677 (2001) [arXiv:astro-ph/0103134]. 
L. Verde et al. [WMAP Collaboration], Astrophys. J. 
Suppl. 148, 195 (2003) [arXiv:astro-ph/0302218]. 
D. MacKay, Information Theory, Inference, and Learn- 
ing Algorithms, Cambridge University Press, Cambride, 
U.K. (2003) 

M. B. Hofltman and M. S. Turner, Phys. Rev. D 64, 
023506 (2001) [arXiv:astro-ph/0006321]. 
A. R. Liddle, Phys. Rev. D 68, 103504 (2003) 
[arXiv:astro-ph/0307286] . 

R. Easther and H. Peiris, JCAP 0609, 010 (2006) 
[arXiv:astro-ph/0604214] . 

M. Liguori and A. Riotto, arXiv:0808.3255 [astro-ph]. 



